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Abstract. We derive the exact bifurcation diagram of the Duffing oscillator with parametric noise thanks 
to the analytical study of the associated Lyapunov exponent. When the fixed point is unstable for the 
underlying deterministic dynamics, we show that the system undergoes a noise-induced reentrant transition 
in a given range of parameters. The fixed point is stabilised when the amplitude of the noise belongs to 
a well-defined interval. Noisy oscillations are found outside that range, i.e., for both weaker and stronger 
noise. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.10.Gg 
Stochastic analysis methods (Fokker- Planck, Langevin, etc.) - 05. 45. -a Nonlinear dynamics and nonlinear 
dynamical systems 



In a classical calculation, Kapitza (1951) has shown 
that the unstable upright position of an inverted pendu- 
lum can be stabilised if its suspension axis undergoes si- 
nusoidal vibrations of high enough frequency 0. More 
generally, stabilisation can also be obtained with random 
forcing In both cases, analytical derivations of the 

stability limit are based on perturbative approaches, i.e., 
in the limit of small forcing or noise amplitudes. 

In a recent work J51, we studied the Duffing oscilla- 
I tor with parametric white noise when the fixed point of 
. the underlying deterministic equation is stable: a purely 
' noise-induced transition |H] occurs when stochastic forc- 
, ing is strong enough compared to dissipation so as to 'lift' 
' the system away from the absolute minimum of the po- 
I tential well. We showed, using a factorisation argument, 
. that the noise-induced transition occurs precisely when 
' the Lyapunov exponent of the linearised stochastic equa- 
tion changes sign. An analytical calculation of the Lya- 
punov exponent allowed us to deduce for all parameter 
values the bifurcation diagram, that was previously known 
only in the small noise (perturbative) limit. The relation 
between stochastic transitions in a nonlinear Langevin 
equation and the sign of the Lyapunov exponent is in fact 
mathematically rigorous and has been proved under fairly 
general conditions [7]. In [3], we also made the following 
striking observation: close to the bifurcation, the averaged 
observables of the oscillator (energy, amplitude square and 
velocity square), as well as all their non-zero higher-order 
moments, scale linearly with the distance from the thresh- 



old. This multifractal behaviour may provide a generic 
criterion to distinguish noise-induced transitions from or- 
dinary deterministic bifurcations. 

In the present work, we reformulate the stochastic sta- 
bilisation of an unstable fixed point of the underlying de- 
terministic system within the framework of noise-induced 
transitions. We extend the analysis of 5 and derive the 
full phase diagram of an inverted Duffing oscillator. We 
show in particular that a reentrant transition occurs in 
this zero-dimensional system. (Reentrant transitions in- 
duced by noise have been studied in the more complex 
setting of spatially extended systems [5|). Moreover, the 
inverted Duffing oscillator also exhibits a multifractal be- 
haviour close to the transition point. 

The dissipative stochastic system considered here is 
defined by the equation: 



d^x da; /— , dU 
d^+^d^ + ^^W^ = -a^ 



(1) 



where x{t) is the position of the oscillator at time t, 7 the 
dissipation rate and ^{t) denotes a stochastic process. The 
confining, anharmonic potential l/({x) is defined as: 



U{x) 



1 2^1 4 
-HX +-X 



(2) 



where /it is a real parameter. Without noise, the corre- 
sponding deterministic system undergoes a forward pitch- 
fork bifurcation when the origin becomes unstable as fj. 
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changes sign from negative to positive values. In 0, we 
only considered the case /i < and showed that for strong 
enough noise, the origin becomes unstable. Here, we study 
the case fi > ('inverted' Duffing oscillator) and show 
that for a finite range of positive values of /i, a reentrant 
transition is observed when the noise amplitude is varied: 
the noisy oscillations obtained for weak and strong noise 
are suppressed for noise of intermediate amplitude. Our 
results are non-perturbative: they are based on an exact 
calculation of the Lyapunov exponent, performed for arbi- 
trary parameter values, ^(t) being a Gaussian white noise 
process. Furthermore, our numerical simulations indicate 
that the phenomenology described above is unchanged for 
coloured Ornstein-Uhlenbeck noise. 

We first rescale the time variable by taking the dissi- 
pative scale 7"^ as the new time unit. Equation then 
becomes 
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Fig. 1. Lyapunov exponent /l of a hnear, damped oscillator 
with parametric white noise. For different values of the control 
parameter a, we plot A vs. the noise amplitude A. Inset: Graph 
of the function g{x) (equation 11311 '). 



d X dx 



■ax + x-^ + ^ A i{t) X = 0, 



dt2 dt 

where we have defined the dimensionless parameters 
a = — and A — — ^ 



(3) 



(4) 



and rescaled the amplitude x{t) by a factor 7. Linearis- 
ing equation (PJ about the origin, we obtain the following 
stochastic differential equation: 

5 + ^-«a; + vze(t)x = o. (5) 

The Lyapunov exponent of a stochastic dynamical system 
is generally defined as the long-time average of the local 
divergence rate from a given orbit [H]. In the case discussed 
here, deviations from the (trivial) orbit defined by the 
origin in phase space, {x{t),x(t)) = (0, 0), satisfy equation 
(jSJ). In practice, we use the (equivalent) definition for the 
(maximal) Lyapunov exponent A 



A = lim 

t — >oc 



(log" 



(6) 



where the brackets denote ensemble averaging. Let z{t) — 
x{t)/x{t). From equation ^ we find that the new variable 
z{t) obeys: 



The stationary p.d.f. of z can then be computed exactly, 
by adapting a calculation presented in detail in |2] (see also 
|l()llllj V where the Fokker-Planck equation equivalent to 
equations l|7|-l|njl is solved in the stationary regime. We 
find: 



Psta.t{z) — — 



exp<j -(if(z) 



^(.y))\dy, (10) 



where tf'(y) = ay — — ^y^, and TV is a normalisa- 
tion constant. After some calculations along the lines of 
Ref. inj, equation © leads to the following expression for 
the Lyapunov exponent: 



yl(a,Z\) 



/c 



+00 dtt 
V^i 



(11) 



We now examine the properties of yl as a function of a 
and A. We shall limit the discussion to the range a > — | 
(see 121 for a detailed study of the range (— cxd,0]). Let 
A = A/{a+ D'^/^. We rewrite equation l(TT1l as: 



1 



1 



Aia, A) = - g{A) + 1^ {VT 



Aa-1) 

where the function g{x) is defined for a; > as 



(12) 



z = a-z-z^-\/A ^{t) . 
The Lyapunov exponent A is equal to 

A = (2;)stat = / -Z Psta,t{z) dz 



(7) 



(8) 



where Psta,t{z) is the stationary probability distribution 
function (p.d.f.) of the variable z, solution of equation Q. 

External fiuctuations acting upon the oscillator will 
first be modeled by Gaussian white noise ^ (t) of zero mean 
value {£,{t)) — and of unit amplitude: 



mm)^5{t~t'). 



(9) 



The method of steepest descent yields 



5(0) = , 



9'iO) ^ 4 



(13) 



(14) 

The function g{x) decreases in the interval [0,x„i], with 
an absolute minimum g{xm) — —0.54 at Xm — 3.40, then 
g{x) increases to infinity over [xrm+oo) (see the inset of 
Fig-D- At fixed a, the behaviour of A{a, A) with respect 
to A is deduced from that of g{x) by translations and 
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Fig. 2. Bifurcation diagram of the Duffing oscillator with para- 
metric white noise. The solid line is the locus in parameter 
space {a, A) where A{a, A) = 0. The bifurcation line a = of 
the noiseless dynamical system is drawn for comparison (dot- 
ted line). 



dilations along coordinate axes. A few representative ex- 
amples are drawn in Fig. ^ (z) a = —0.2: A changes sign 
once (noise-induced bifurcation of the nonlinear system 
0); (m) a — 0.2: A changes sign twice (noise-induced reen- 
trant transition of the nonlinear system); (iii) a = 0.5: A 
is positive for all A (no transition). 

From this analysis, the bifurcation diagram of the non- 
linear oscillator with parametric, Gaussian white noise 
(equations (jniandEJ) follows immediately. The origin is 
stable (resp. unstable) when the Lyapunov exponent of 
the linearised system is negative (resp. positive). The bi- 
furcation line a — ac{A) is defined by the equation 



Aia,{A),A)^0. 



(15) 



The transition is best qualified as a stochastic Hopf bi- 
furcation: the bifurcated state displays noisy oscillations 
around the origin, with a non-zero r.m.s. of the oscilla- 
tor's position and velocity. For A < A* ~ 3.55 (resp. 
A > A*), the origin is stabilised (resp. destabilised) by 
the stochastic forcing in the range < a < ac{A) (resp. 
ctc{A) < a < 0). The full bifurcation diagram of the in- 
verted stochastic Duffing oscillator is displayed in Fig. |21 
In the weak noise limit Z\ ^ 0+, we obtain from Eq. H14|l 
ac{A) ~ (-2 g'{0) A) ^ A/2, in agreement with the 
result of the (perturbative) Poincare-Linstedt expansion 
performed in 

We emphasise that the stability of the origin of the 
nonlinear random dynamical system Q cannot be de- 
duced from a stability analysis of finite-order moments of 
the linearised system |12) : indeed second-order moments 
of solutions of equation |(SJ with white noise forcing are 
always unstable when a is positive. The proper indicator 
of the transition of the nonlinear system is the Lyapunov 
exponent of solutions of the linearised equations. 

As an example, we show in Fig. |2| numerical evidence 
of a reentrant transition observed for a — 0.2. The asymp- 
totically stable state is the origin when the noise ampli- 
tude belongs to the bounded interval A G [Ai, A2], where 
Ai and A2 are the two solutions of A{a = 0.2, A) = 0. 
Noisy oscillations are found for both weaker (A < Ai) 




Fig. 3. Reentrant transition of the noisy nonlinear oscilla- 
tor ^ subject to parametric white noise lO, with a fixed 
parameter a — 0.2. The numerical solutions of the equation 
yl(0.2,Zi) = are Ai ~ 0.627, A2 ~ 1.700. Even-order mo- 
ments (a;^), (a;'*) and (x^) are plotted versus the distance to 
threshold e — {Ai — A)/Ai (symbols). Dashed lines respecting 
a linear behaviour (a;^") cx are drawn to guide the eye. Lin- 
ear scaling is also observed for the second transition {A > A% 
not shown). Inset: the mean square position (measured in the 
stationary regime) is non-zero for A £ [0, Zii[U]zi2, +oo[. 



and stronger {A > A2) noise strengths. As predicted in 
|S] , even-order moments of the position and velocity scale 
linearly with the distance to threshold in the vicinity of 
these stochastic bifurcations. (Note that odd-order mo- 
ment are equal to zero by symmetry.) 

Since the locus of the transition is determined by a 
dynamical property of the linearised system, it cannot de- 
pend on the precise functional form of the confining po- 
tential L{{x) for large x. For instance, the bifurcation line 
is unchanged when the confining potential reads: 



U{x) = —^ax 



6 



(16) 



i.e., when the deterministic nonlinear system undergoes 
a backward pitchfork bifurcation. We checked numerically 
that such is indeed the case. 

Let us now turn to stochastic forcing by coloured noise. 
For definiteness, we shall use an Ornstein-Uhlenbeck pro- 
cess £,{t) of correlation time r, defined as the solution of 
the stochastic differential equation: 



(17) 



where r]{t) denotes Gaussian white noise with zero mean 
and unit amplitude [{'r]{t)ri{t')) = 5{t — t')). The analytic 
expression of the stationary (marginal) p.d.f. Pstat(-z) as- 
sociated with the linear system ()7ll7ll is not known: the 
Lyapunov exponent must be evaluated numerically. The 
line a = ac{A) with A{ac{A), Z\) = is drawn in Fig.^a 
for r = 0.1 and 1.0. (We checked that the same measure- 
ment protocol, when used with white noise forcing, yields 
data in agreement with the analytic result Hll|l .) Numer- 
ical simulations confirm that, with parametric coloured 
noise, (z) bifurcations of the nonlinear system (equations 
(|3ll7|l ) occur where the Lyapunov exponent of the linear 



4 



Kirone Mallick, Philippe Marcq: Noise-induced reentrant transition of the stochastic DufBng oscillator 



(a) 





' 1 


1 1 1 1 1 1 
A > 










A < 






— T = 0.0 






-- T = 0.1 






■-■ T = 1.0 






1 


I.I.I, 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 



A 
(b) 



A 0.05 



-0.0: 





<:^'z = 0.0 
Q-e-u = 0.0 

Ur^Z = 1.0 
B-ET = 1.0 


• T^T-^ ^ — 7" 

/ p / / 
r/ / ^ 


p' 



j^L^^ \ — \ , \ , 1 

V.O 0.1 0.2 0.3 0.4 



a 

Fig. 4. Parametric forcing with Ornstein-Uhlenbeck noise of 
correlation time r. (a) Bifurcation lines for r — 0.1 and 1.0 
are obtained from numerical measurements of the Lyapunov 
exponent @. For comparison, we draw the (analytic) white 
noise line already given in Fig. El (6) For r = 1.0 and 0.0 
(white noise), we plot on the same graph, for a fixed value 
A = 1.096, and with respect to a, the numerical values of the 
Lyapunov exponent of the linear system (symbols with dashed 
lines) and of the averaged square position in the oscillatory 
regime of the nonlinear system (symbols with solid lines) . The 
bifurcation occurs where the Lyapunov exponent changes sign. 
Note the linear dependence of {x^) on distance to threshold. 



system changes sign; (ii) averaged observables scale lin- 
early with distance to threshold close to the bifurcation 
(see Fig. Ejb). In the weak noise limit, our data agrees 
with the prediction of E]: a^A) - A/{2{1 + t)). For a 
noise amplitude of order 1, the bifurcation hnc is quaUta- 
tively similar to that obtained with white noise. However, 
depending on the value of A, the value of the bifurcation 
point ac{A) is not necessarily a monotonic function of t. 

We have shown here that noise can suppress oscil- 
lations by stabilizing a deterministically unstable fixed 
point. Our study can be related to similar observations 
in the context of chaotic systems where Lyapunov expo- 
nents also play a crucial role as indicators of transitions 
|13j . Indeed, the effective lowering of Lyapunov exponents 
by noise can induce synchronization in a pair of chaotic 
systems P^ItT^ . or suppress chaos in nonlinear oscillators 
[W\ or in neural networks |17j. 



Our hope is that this work will help to bridge the gap 
between the abstract theory of random dynamical systems 
[7| and experimental investigations of the stability of phys- 
ical systems subject to random forcing |18II19| . The be- 
haviour of Lyapunov exponents of the linearised stochastic 
system (or "sample stability" ) is known to explain the sta- 
bility properties of a regular state of electrohydrodynamic 
convection of nematic liquid crystals driven by multiplica- 
tive, dichotomous noise |20| . A non-perturbative study of 
bifurcations of spatially-extended systems under stochas- 
tic forcing of arbitrary strength |21j remains a challenging 
open problem. 



References 

1. L.D. Landau, E.V. Lifchitz, Mechanics (Pergamon Press, 
Oxford, 1969) 

2. R. Graham, A. Schenzle, Phys. Rev. A 26, 1676 (1982) 

3. M. Liicke, F. Schank, Phys. Rev. Lett. 54, 1465 (1985); M. 
Liicke, in Noise in Dynamical Systems, Vol. 2: Theory of 
Noise-Induced Processes in Special Applications, edited by 
F. Moss, P.V.E. Mc Clintock (Cambridge University Press, 
Cambridge, 1989) 

4. P.S. Landa, A.A. Zaikin, Phys. Rev. E 54, 3535 (1996) 

5. K. Mallick, P. Marcq, Eur. Phys. J. B 36, 119 (2003) 

6. H. Horsthemke, R. Lefever, Noise Induced Transitions 
(Springer- Verlag, Berlin, 1984) 

7. L. Arnold, Random Dynamical Systems (Springer- Verlag, 
Berlin, 1998) 

8. C. van den Broeck, J.M.R. Parrondo and R. Toral, Phys. 
Rev. Lett 73, 3395 (1994) 

9. L. Schimansky-Geier, H. Herzel, J. Stat. Phys. 70, 141 
(1993) 

10. L. Tessieri, F.M. Izrailev, Phys. Rev. E 62, 3090 (2000) 

11. P. Imkeller, C. Lederer, Dyn. Systems 16, 29 (2001) 

12. B.J. West, K. Lindenberg, V. Seshadri, Physica A 102, 470 
(1980); K. Lindenberg, V. Seshadri, B.J. West, Phys. Rev A 
22, 2171 (1980); Physica A 105, 445 (1981) 

13. A.S. Pikovsky, Phys. Lett. A 165, 33 (1992) 

14. A. Maritan, J.Y. Banavar, Phys. Rev. Lett. 72, 1451 
(1994); A.S. Pikovsky, Phys. Rev. Lett. 73, 2931 (1994) 

15. E. Sanchez, M.A. Mafias, V. Perez-Munuzuri, Phys. Rev. 
E 56, 4068 (1997) 

16. S. Rajasekhar, M. Lakshmanan, Physica A 167, 793 
(1990); Physica D 67, 282 (1993) 

17. L. Molgedey, J. Schuchhardt, H. G. Schuster, Phys. Rev. 
Lett. 69, 3717 (1992) 

18. R. Berthet, S. Residori, B. Roman, S. Fauve, Phys. Rev. 
Lett. 33, 557 (2002); R. Berthet, A. Petrossian, S. Residori, 
B. Roman, S. Fauve, Physica D 174, 84 (2003) 

19. F. Petrelis, S. Aumaitre, Eur. Phys. J. B 34 281 (2003) 

20. U. Behn, A. Lange, T. John, Phys. Rev. E 58, 2047 (1998) 

21. A. Becker, L. Kramer, Phys. Rev. Lett. 73, 955 (1994); 
Physica D 90, 408 (1996); J. Roder, H. Roder, L. Kramer, 
Phys. Rev. E 55, 7068 (1997) 



